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AUTOMATIC GUIDANCE AND CONTROL OF A TRANSPORT AIRCRAFT 
DURING A HELICAL LANDING APPROACH 

Daniel J. Crawford and Roland L. Bowles 
Langley Research Center 

SUMMARY 

Linear optimal regulator theory is applied to a nonlinear simulation of a transport 
aircraft performing a helical landing approach. A closed-form expression for the quasi - 
steady nominal flight path is presented along with the method for determining the corre- 
sponding constant nominal control inputs. The Jacobian matrices and the weighting 
matrices in the cost functional are time varying. A method of solving for the optimal 
feedback gains is reviewed. The control system is tested on several alternative landing 
approaches using both 3° and 6° flight path angles. On each landing approach, the air- 
craft was subjected to large random initial-state errors and to randomly directed cross- 
winds. The system was also tested for sensitivity to changes in the parameters of the 
aircraft and of the atmosphere. Performance of the optimal controller on all the 3° 
approaches was very good. Mean errors in the lateral, vertical, and longitudinal direc- 
tions for the principal flight path (270° turn) were 0.8 m, -0.1 m, and -9.4 m. The magni- 
tude of the respective ranges of these errors were 1.0 m, 1.9 m, and 26.7 m. The con- 
trol system proved to be reasonably insensitive to parametric uncertainties. Perform- 
ance was not as good on the 6° approaches. On the principal flight path, mean errors for 
the lateral, vertical, and the longitudinal directions were 9.4 m, -0.4 m, and -2.9 m and 
the magnitude of the respective ranges for these errors were 8.7 m, 3.2 m, and 45.2 m. 

A modification to the 6° flight path is proposed for the purpose of improving performance. 

INTRODUCTION 

Aircraft traffic in the neighborhood of commercial airports has been the subject of 
intensive study during the past several years. Noise, collision avoidance, airspace con- 
gestion, air pollution, and dangerous wing -tip vortices are all problems which are being 
studied by engineers today. Presently, in the near-terminal area, aircraft typically fol- 
low essentially straight-line paths to a common point, where they are sequenced to follow 
one another down a 3° flight path on a straight-in approach to the runway. This procedure 
very often results in low-altitude flying at relatively high power settings over residential 
or densely populated business districts which is objectionable, at least, from a noise 


standpoint if not from the standpoint of safety. Recently, steeper, 6° flight paths and 
two-segment final approaches (ref. 1) have been proposed to eliminate these problems 
and are currently being tested. This paper proposes an alternative to these procedures. 

It consists of descent along a helical path to a low altitude where the straight-in final 
approach is intercepted and followed. This flight path has the advantage of keeping the 
aircraft at low power and at relatively high altitude except in the immediate area of the 
airport. The approach path is helical, the flight -path angle, ground speed, and turning 
rate are maintained constant during descent, and the control inputs are fixed except to 
null flight -path errors. 

Implementation of the proposed automatic guidance system requires *wo factors 
which are not at present universally available. The two factors are an on-board flight 
control computer and a fairly precise knowledge of the aircraft's position in space. 
However, with the advent of automatic landing systems, flight computers are becoming 
much more common on commercial transports and the proposed installation of the micro- 
wave landing system (MLS) equipment at airports throughout the country will make the 
necessary position data available (ref. 1). 

The method proposed in this paper uses the theory of optimum control which is 
treated in references 2 to 6. The theory is applied to a three -degree -of -freedom non- 
linear aircraft simulation. In order to pose the problem as a linear regulator, the equa- 
tions were linearized about a nominal state (helical path) and control trajectory. The 
resulting near -optimum feedback control gains were then tested in the nonlinear model. 
This work is an extension of that presented in references 7 and 8. The extension includes 
increasing the number of states to six and the number of controls to three. The additional 
states remove the restriction that altitude be monotonically decreasing and the inclusion 
of thrust (autothrottle) as a control allows treatment of aircraft other than gliders. Ref- 
erence 7 concerns itself with a HL-10 lifting body reentry vehicle and reference 8 with an 
unpowered Boeing 707 with no external wind disturbances. Another extension of this work 
is inclusion of the quasi-steady helical path with its consequent piecewise constant control 
inputs. This factor greatly reduces the problem of computer storage of the state and con- 
trol trajectories. The method for computing these quasi -steady paths along with the 
method for computing appropriate control inputs are presented herein. In addition, a 
derivation of the equations of motion is given, along with the method of solving for the 
near -optimum feedback gains. It should be noted that this study, as do those of refer- 
ences 6 and 8, uses angle of attack and bank angle as control inputs. This is a simplifi- 
cation of the dynamic model which does not include control surfaces such as elevators, 
ailerons, or rudder. 

Feedback gains were obtained for several candidate approach paths including both 
3° and 6° flight -path angles and turns of 180°, 270°, and 360°. An overhead approach 
was the principal path studied (i.e., one in which the aircraft crossed over the runway at 
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a 90° angle heading east at an altitude of 550 m). At about 1800 m east of the runway 
centerline it intercepted the helix and descended along it, making a total turn, to the right, 
of 270°. The turn was completed at an altitude of 100 m where the aircraft intercepted 
the straight -in approach. At this time, it was about 1900 m from the touchdown point on 
the runway. The plane continued to descend until at a 12-m altitude, the onset of the 
flare maneuver, the simulation was terminated. The effectiveness of the control system 
was tested in the nonlinear simulation by introducing initial errors in the states, cross- 
winds, and uncertainties in the aerodynamics of the aircraft. 

SYMBOLS 

A, A(t) state -coefficient matrix (6 x 6) in linearized equations of motion 

B, B(t) control-coefficient matrix (6 x 3) in linearized equations of motion 

Cp> total drag coefficient 

Cp) 0 drag coefficient for zero lift 

Cp j total lift coefficient 

c L a = aC L l ba > rad_1 

D aerodynamic drag force, N 

f ,F(x,u) derivative of the state vector x (6x1) 

fp vector function (3 x 1) used in determining nominal open-loop control for 

equilibrium 

g gravitational force per unit mass, 9.80 m/sec^ 

h altitude (-z), m 

H Hamiltonian of linearized system 

J cost functional (performance index) 
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K,K(t) 


L 


m 

M 


P,p(t) 

q 

Q 

r 

R 

S 


Si,S!(t) 


t 

T 


T/m 

Ty 


u 


feedback -gain matrix (3 x 6) 
i,j element of K 

aerodynamic lift force, N 
mass of aircraft, kg 

weighting matrix (6 x 6) on terminal-state variations 

vector (6 x 1) of costates 

dynamic pressure, N/m^ 

weighting matrix (6 x 6) on state variations 

radius of helix, m 

weighting matrix (3 x 3) on control variations 
reference wing area of aircraft, m^ 

matrix (6 x 6) relating the costates to the states in equation (26), solution of 
the time -varying matrix Riccati equation 

time, independent variable, sec 

total engine thrust. N 

thrust to mass ratio, m/sec^ 

orthogonal transformation (3 x 3) through the angle y 
orthogonal transformation (3x3) through the angle 0 
orthogonal transformation (3 x 3) through the angle p 
control vector (3 x 1), [a 0 T/m]^ 
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V 


magnitude of the velocity vector with respect to the inertial axis system, 
m/sec 

V velocity vector (3 X 1) resolved along rotating axes system 

V w scalar magnitude of crosswinds, m/sec or knots 


x 


x 

X, Y,Z 


X’.Y’.Z’ 


y 


z 


z 

Z,Z(t ) 
a 


a o 


displacement of aircraft in the north-south direction relative to the origin of 
the fixed axes system, m 

state vector (6 x 1) of the system, [x y z V y it']'*’ 
inertial axis system 
rotating axis system 

displacement of aircraft in the east -west direction relative to the origin of 
the fixed axes system, m 

vertical displacement of aircraft relative to the origin of the fixed axes 
system, m 

sink rate, m/sec 

state coefficient matrix (12 x 12) in equation (23b) 

angle of attack, deg 

angle of attack for zero lift, deg 


y flight -path angle, rad or deg 

6u variation from nominal control vector (3 x 1) 

6x variation from nominal state vector (6 x 1) 


8q!r, 50 r , 

S(T/m) R 


symbols used in defining the diagonal elements of the R weighting matrix 
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6 x M» 6 YM^ 

8 z m , 8 Vm \ symbols used in defining the diagonal elements of the M weighting matrix 

^M.^M 

J 

6 x Q’ 6 Yq" 

^ z Q>^Vq v symbols used in defining the diagonal elements of the Q weighting matrix 

6y Q’ 6 ^Q 

J 

77 aerodynamic efficiency factor in equations (6) 

p,p(h) density of air, kg/m 3 

<fi aircraft roll angle, rad or deg 


if/ aircraft track angle, rad or deg 

if/ w direction of wind, rad or deg 


u> 

tt,fi(t,h) 



rotation vector of aircraft (3 x 1 ) resolved along rotating axis system 

transition matrix (12 x 12) for the combined system of states and costates in 
equation ( 24 ) 

submatrices (6 x 6) of fi 


Subscripts: 

0 initial value 

f final value 

N nominal value 

SI nominal value when the aircraft executes roll-out maneuver 
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Operators: 


a x b cross product of a and b 

U[a,b] uniformly distributed random variable with range (a,b) 

[ ]“* inverse of matrix 

[ ] T transpose of matrix or vector 

A dot over a variable denotes the derivative with respect to time. 

ANALYSIS 

The purpose of this section is to present the mathematics involved in simulating the 
aircraft and in implementing the feedback control law. As illustrated in figure 1, the air- 
craft is modeled by a set of six nonlinear ordinary differential equations represented by 
the single vector equation 

x = F(x,u) (1) 

This set of equations was integrated numerically using a fixed-step fourth-order Runge 
Kutta algorithm. The feedback control law used in this study requires that the nominal 
state trajectories x^, the nominal control time histories u^, and the feedback gain 
matrix K be stored in an on-board computer. However, it will be shown that by choos- 
ing particular nominal trajectories, the vector u^ changes only twice, and the vector 
x N is a simply computed analytic vector function of time. This eliminates some of the 
on-board storage requirements and makes the control system simpler and more flexible. 

As simulated, the aircraft has six state variables and three control inputs defined 
as 


x = [x y z V y xp]^ 
u = [a <fi T/m] T 

The state vector x is made up of the variables x, y, and z describing the aircrafts 
position; V, its inertial speed; y, its flight -path angle; and ip, its track angle. The 
control vector u is made up of the angle of attack a?, the roll angle 0, and the ratio of 
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engine thrust to aircraft mass T/m. The two perturbation vectors 6x and 6u, shown 
in figure 1, are defined as 

<5x = x - x-^ 

5u = u - u-^ 

However, in practice, Su is computed as: 

6u = KSi (2) 

where K is a 3 x 6 time -varying gain matrix. This gain matrix is the solution of the 
regulator problem which results from linearizing equation (1) about the nominal flight 
path and control vector, and it is precomputed for a particular nominal flight path. This 
linearization yields: 

6x = A(t)6x + B(t)6u = [A(t) + B(t)K(t)]6x (3) 

where A(t) = 9f /9x evaluated at x = x^ and u = ujj, B(t) = 9f /9u evaluated at 
x = x N and u = u N , A(t) is a 6 x 6 matrix, and B(t) is a 6 x 3 matrix. A statement 
of the regulator problem is then to find K in equations (2) and (3) which minimizes 
deviations of the states and controls from their nominal values. 

More formally, the problem is to determine K which minimizes the quadratic 
performance index J given by 


J = \ 5x T (t f )M6x(t f ) 



6x(t)^Q6x(t) + 5u(t)^R<5u(t) 


dt 


( 4 ) 


subject to the constraints of equation (3) and where M, Q, and R are weighting 
matrices. 


Equations of Motion 

The inertial axis system illustrated in figure 2 has its origin at the glide -path 
intercept point on the runway. The X-axis is parallel to the runway and is positive in 
the direction of landing. The Y-axis is horizontal, perpendicular to the X-axis and is 
positive to the right. The vertical Z-axis is perpendicular to both X and Y, and is 
positive downward. The aircraft is treated as a point mass, subject to the forces of 
gravity, engine thrust, and aerodynamic lift and drag. The rotating axis system associ- 
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ated with the aircraft is chosen such that its X' -axis is always alined with the aircraft's 
velocity vector. Two Euler angle transformations relate the inertial axis system to the 
rotating axis system. They are 


T 




COS 1 p 
-sin xp 


0 


sin p 0 

cos p 0 

0 1 


and 


cos y 


T y = 


o 


-sin y 


0 sin y 

1 0 

0 cos y 


These rotations are shown in figure 2. This axis system does not roll with the aircraft. 
Therefore, the Y' -axis remains horizontal and the Z'-axis stays in a vertical plane with 
the velocity vector. The rotation vector resolved along the rotating axis system is 

To =[p sin y -y p cos y]"^ 


and the velocity vector is 
V = [V 0 0] T 

The forces acting on the aircraft resolved along the rotating axes are 

Drag = [-D 0 0]"^ 

Thrust = T^[T cos a 0 -T sin of]'*’ 

Gravity = TyTj^O 0 rag]^ 

Lift = T^[0 0 -L] t 
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where 


T 


0 = 


1 

0 

0 


0 

cos 0 
sin 0 


0 

-sin 0 
cos 0 


The lift force is perpendicular to the velocity and rolls with the aircraft. The transfor- 
mation T ^ resolves this force along the Y'- and Z' -rotating axes. The angle of attack 
a is defined as the angle between the velocity vector and the thrust vector of the air- 
craft. The dynamics of the aircraft are then described by equating the change in linear 
momentum to the sum of the applied forces 

m(V + cZ> x V) = Drag + Thrust + Gravity + Lift 


or 



V 


0 


’-d" 


T cos a 


mg sin y 


0 

m 

0 

+ m 

V0 cos y 

= 

0 

+ 

T sin a sin 0 

■f 

0 

4- 

L sin 0 


0_ 


Vy 


_ 0 _ 


-T sin a cos 0 


mg cos y 


-L cos 0 


To facilitate integration of these equations, they are written with the derivatives of the 
states V, 0, and y on the left-hand side. These three equations are combined with 
three kinematic relations which resolve the velocity along the X-, Y-, and Z -inertial 
axes 


x 


z 



-1 


V 

0 

0 


The six equations of state are then 


x = V cos y cos 0 

(5a) 

y = V cos y sin 0 

(5b) 

z = V sin y 

(5c) 
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(5d) 


• "P\ »X> 

V = — + — cos a + g sin y 
mm 


y = ^ cos y 


cos 0 - 

mV 


mV 


sin a cos 0 


(5e) 


^ _ (L + T sin a) sin 0 
v mV cos y 


(5f) 


Equations (5) are the expanded form of equation (1). 

Aerodynamics and Aircraft Parameters 

The aircraft simulated in this study is a small two-engined transport airplane 
designed primarily to operate from short runways over relatively short distances. It is 
described in references 9 and 10. 

The equations for the lift and drag forces are 
L = qSC L 


D = qSC D 


where S is the wing reference area and C L and C D are the coefficients of lift and 
drag, respectively, and q is dynamic pressure which are defined as 


C L = c L a (“ " «o) 


180 


C D = c D,o + 7 ? C L / 



( 6 ) 


where a parabolic drag polar has been assumed. This assumption is not restrictive for 
the range of speed and angle of attack employed in this study. Air density p was 
approximated by the following equation taken from reference 8. 


p = 1.22 


1 - (2.257 x 10" 5 )h 


4.255 


(7) 
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This proved to be a very good approximation to the standard atmospheric density of ref- 
erence 11. The airplane was configured with the landing gear down and the flaps 
deployed at 40°. The parameters of equations (6) a °> ^D,o> anc * V) were 

chosen to fit data on Cl and provided by the manufacturer for this configuration. 

The parameters of the airplane are listed in table 1. 

Nominal Flight Paths 

Linear regulator theory does not in any way restrict the choice of nominal flight 
path. However, practical considerations including ease of implementation do suggest 
certain types of landing approaches. If a path on which the derivatives of velocity and 
flight -path angle are identically zero is chosen, then the nominal state trajectories can 
be expressed as simple analytic functions of time. In which case, much less on-board 
computer memory is needed to store the path. Before examining the helical landing 
approach, first look at the conventional 3°, nondecelerating, straight-in approach. A 
mathematical description of this path is 

x = [V 0 cos y 0 0 V Q sin y Q 

x(0) = [x 0 0 z 0 V 0 y 0 

As the aircraft approaches the ground, it leaves this quasi-steady condition and flares in 
order to touch down smoothly. This path calls for constant control inputs varying only 
to offset atmospheric disturbances. Because of its simplicity, it can be easily 
implemented. 

In order to relieve airport noise problems and to decrease traffic congestion in the 
neighborhood of the airport, various other approaches have been proposed and are being 
tested. These include the steeper 6°, nondecelerating, straight-in approach and a com- 
bination of the 6° and 3° approach. The method proposed in this paper is a constant- 
control helical approach. All other constant -control paths are special cases of this gen- 
eralization. That is, the nominal path is descent along a helix from some initial altitude. 
The centerline of the runway is tangent to the helix, and at some specified lower altitude, 
the aircraft follows this tangent into a straight-in approach. By making the assumption 
that atmospheric density is constant over the range of altitude considered, this path calls 
for fixed controls during the helical descent, and fixed (but different) control inputs dur- 
ing the straight-in portion. The variation in atmospheric density in this range of altitude 
(0 to 700 m) is about 4 percent. The assumption of constant density is made during cal- 
culations of nominal path, nominal control inputs, and feedback gains. During the simu- 
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lated flights, using the closed-loop control system, atmospheric density is allowed to 
vary according to equation (7). The small error introduced by this assumption was 
nulled out by the control system as though it were an external disturbance. A more 
formal description of the nominal path is 


Vq COS Yq COS <// 


V 

Vq cos Yq sin 4? 


y 0 

Vq sin y Q 

x(0) = 

z 0 

0 


v o 

0 


y 0 

1 

o 


_*o 


for the altitude range 

- z o = h = - z si 

where -zgj is the altitude at which the aircraft rolls into the straight -in approach. The 
parameters zq, zgi, Vo, Yq, and i^q are chosen with some degree of latitude and 
they essentially determine the nominal flight path. The last three parameters (x Q , y Q , 
and i//q) are functions of the former and are chosen such that x, y, and xjs have the 
desired values when z equals Zgj. The subscript SI will be used to specify the 
value of a state when z equals zgj. Because the origin of the axis system is at the 
glide-path intercept point on the runway, the desired values of \{/ gj and ygj are zero. 
In order to determine x gI , it is necessary to look at the nominal path after it starts the 
straight -in approach. Its description is 

x = JVq cos 7 q 0 Vq sin Yq 0 

x(t S l) = [ X SI 0 Z SI V 0 Y 0 
During the straight-in part of the flight x and 

<k = iZi=coty 0 

dz dz/dt 0 



z are constant and consequently, 
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is constant. In order to touch down at the origin, xgj and zgj must be related by 


or 


x si = z si f 


X SI = 


Z SI 

tan y Q 


(10) 


Also, during the helical descent, the rate of turn ^ and the rate of descent z are con- 
stant. In order to have the aircraft flying in the correct direction (\p = 2n) when z 
equals zgj, it is necessary that 


d^ _ 2?r - 
dz Z SI - z 0 

and since 

di^ _ di/'' dz 
dt dz dt 

it follows that 

. /277 - tf/ 0 \ 

ipn = ( ^ V 0 sin y n 

0 \ Z SI - z 0/ ° 0 


(ID 


( 12 ) 


The axis of the helix is vertical with origin at 


(x,y) = ( x si> r ) 


where the radius r is given by 


r = 


tan Yq ^ 


d^ 

dz 


(13) 


and through geometrical considerations 


x 0 = x SI + r sin 


(14) 


14 



y 0 = r (i - cos p Q ) 


(15) 


A short example may clarify these computations (see fig. 3, flight path 2). Let the land- 
ing direction be north = 2 ir) and the aircraft's original heading (t = 0) be east 
(i^O = 7r /2). The velocity is 62 m/sec and the flight -path angle is 6° or 67T/180 rad. 
Initial altitude is 550 m. The aircraft descends along a helical path to a 100 -m altitude 
and then down a straight path to the runway. In this example, 

zq = -550 m 
Z SI = -100 m 
V 0 = 62 m/sec 


y o = Tk 




X 0 = Xgj + r sin q ~ -951 + 909 sin ^ ~ -43 m 
y 0 = r^l - cos = r ~ 909 m 
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Integration of equations (8) gives the following state trajectories along the helix. 



(16) 


Likewise, integration of equations (9) yields the state trajectories along the straight por- 
tion of the flight path 


x = 


x si 


Vo cos Yq 

0 


0 

Z SI 

+ 

v 0 sin ^0 

v o 


0 

^0 


0 

0 


0 


(* - tsi) 


(17) 


The time for the wings -level maneuver tgj is given by 


= Z gI -_Zq_ 

SI V 0 sin y Q 


(18) 


In summary, for a helical approach there is no need to integrate and store the equa- 
tions of motion. Given a particular input parameter set (z 0 , z SI , Vq, Yq, ^q) equa- 
tions (16) and (17) yield the value of the states at any time. Five secondary parameters 
( X SI> &0’ x 0> yo> tsi) necessary to evaluate equations (16) and (17) are given as functions 
of the input parameters by equations (10) to (15) and equation (18). 
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Nominal Control Inputs 


p 


After specifying this particular class of nominal flight paths, it is necessary to 
solve for the appropriate open-loop control. Since, as has been stated earlier in equa- 
tions (8) and (9), 

V = Y = 0 

* = for h S -z gI 

ij/y — 0 for h = ~Zgj 

u is a constant vector on the spiral and a second constant on the straight -in portion. In 
order to make equations (5) agree with equations (8), the following equations must be 
satisfied: 

D N T • A 

— — - — cos a - g sin y n = 0 
mm & u 

T - - sin a cos <p - rh cos Vn + — 77 - cos 4> = 0 
mV 0 V 0 u mV 0 

(Lm + T sin a) sin <p 
rp . = 0 

u mVo cos Vo 

In these equations, 

L N = q N sc L 
d N = qN SC D 

where 

q N = | ^n v o 2 



17 



and Cl and Cj) are defined in equations (6). Using the assumption of constant air 
density discussed earlier, an intermediate value of altitude hjj is chosen, and is 

computed using equation (7). The open-loop control 

u = [or <p T/m]*^ 

which is the solution of equations (19), is not unique and these equations have no closed- 
form solution. Equations (19) are of the form 

fl(u) = 0 


and it was solved by the Newton -Raphs on iterative method described in reference 12. 
This method usually requires the Jacobian matrix which is given by 


9D N 

ao? 


+ T sin a 
m 


0 


-cos a 



T cos a + 


aL N \ 

da ) 


cos (p 


mV 0 


T cos a + sin (p 
mVQ cos Yq 


-(L n + T sin a) sin <p cos <p sin a 
mV 0 V 0 

-(Lr + T sin aj cos (p _ sin a sin (p 
mV 0 cos Yq Vq cos y Q 


where 


8L N _ n e x p T 
acr q N b 180 L ' L Q! 

and 

The possibility always exists that this algorithm will not converge at all or will 
converge to an undesirable value of u. Therefore, the solution must be examined sub- 
jectively before being accepted. A practical way to avoid this problem is to start the 
iterative procedure close to the desired final solution. A suggested method of doing this 
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I 


is to make the assumption that engine thrust is alined with the velocity vector of the air- 
craft. With this assumption, which is in error by the size of the angle of attack, equa- 
tions (19) become 


Dn 

m 


X . g sin r 0 = 0 


^ cos r 0 


+ 


l n 

mV 0 


COS 0 = 0 


• _ Ln sm_0__ = 0 
o mV 0 cos yp 


These equations have the following closed-form solution which can be used as a starting 
point in the Newton -Raphs on algorithm. The solution is 
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In order to solve for the open -loop controls for the straight-in portion of the nom- 
inal flight path, the following equations must be satisfied: 
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Again, the Newton -Raphson algorithm is used to solve the two equations for a and 
T/m. The roll angle 0, of course, is zero. Using the same assumption as before that 
thrust and velocity are codirectional, the following equations can be used to start the 
iterative procedure. 
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Notice that D in equations (20) will have a different numerical value than in 

equations (21). These equations used for starting the Newton -Rap hs on iteration proce- 
dure yield values which are very close to the final values. 

Feedback Control Gains 

The purpose of this section is to outline the solution of the linear regulator prob- 
lem. References 1 to 5 all are concerned, in part, with this problem. A statement 
of the problem is given by equations (1), (3), and (4), which are repeated here for 
convenience. 


x = f (x,u) 


6x = A(t)Sx + B(t)6u = (~A(t) + B(t)K(t)]5x 


J = \ 6x T (tf)M5x(tf) + | ^ f 


5x(t)^Q6x(t) + 6u(t)^R6u(t) 


dt 


A(t) and B(t) are defined as 
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= af_ 
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au 


and M, Q, and R are weighting matrices which must be selected by the control- 
system designer. The problem is to find the matrix K in equation (3) which minimizes 
the quadratic performance index of equation (4). As shown in the references, the state 
equations are combined with the costate equations to form the system: 
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where p is the vector of costates. 1 Equation (23a) can be rewritten in the form 



(23b) 


and if we assume that Z is constant over a small interval of time At then we can 
write the transition equation of the system as 


6x(t + At) 
p(t + At) 
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5x(t) 

p(t) 
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(24) 


where 


S7(t, At) = e Z ^ At = 

fill 
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n 21 
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(25) 


It is further shown that the costates are linearly related to the states; that is, 
p(t) = Sj6x(t) 

and that the final condition on Sj is given by 

Sl(tf) = M 


Combining equations (24), (25), and (26) yields a transition equation for Sj 

= 12(t, At) 


6x(t + At) 

Si(t + At)6x(t + At) 


6x(t) 

Si(t)Mt) 


S x (t + At) [fin + n 12 Si(t)]6x(t) = [n 2 i + ^22^l] 6x(t) 


(26) 


(27) 


iReference 2 shows that M, Q, and R must be symmetric; that M and Q 
must be at least positive semidefinite; and that R must be positive definite. 
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or 


Si(t + At) = [p21 + fi 22 s l(t)] [ fi ll + fi 12 s l(t)] 1 (28) 

Reference 2 proves that this inverse does exist. A comprehensive development of equa- 
tion (28) can be found in chapter 7 of reference 13. 

Using the difference equation (28) with the final value of the matrix Sj given by 
equation (27), we can proceed backward in time and obtain Sj for all t in the range 
(t(j,tf). At each point in t>me, it is necessary to compute 0(t,At) as a function of A(t), 
B(t), R, and Q, using equations (23) to (25). The gain matrix K(t) is obtained by the 
following relationship: 


6u = -R-^SiSxft) 


or 


K(t) = -R-^Si 


(29) 


Implementation 

The landing -approach guidance scheme proposed in this paper would be implemented 
in the following manner: 

(1) For a particular nominal flight path, the control -system designer, after exten- 
sive analysis and testing, chooses weighting matrices M, Q, and R for the perform- 
ance index in equation (4). 

(2) The state trajectories on the nominal flight path are computed using equa- 
tions (10) to (18). 

(3) The nominal open-loop control inputs are computed using the Newton-Raphson 
algorithm to solve equations (19) and (21). This iterative procedure is started with the 
values given in equations (20) and (22), respectively. 

(4) The Jacobian matrices A(t) and B(t) used in equation (3) are computed 
according to equations (Al) and (A2) in the appendix. Equation (3) represents the lineari- 
zation of equation (1) about the nominal trajectory and control. 

(5) The feedback gain matrix K is computed as a function of time using equa- 
tions (23) to (29). 
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(6) In the computer aboard the aircraft are stored the time histories of the nominal 
flight path, the nominal control inputs, and the feedback gains. Alternatively, the algo- 
rithm (eqs. (16) and (17)) which computes the state trajectories is stored. 

(7) The aircraft, or simulated aircraft, is flown into a window (i.e., a region of 
state space which is "close" to the initial states of the nominal flight path). 

(8) When the plane enters the window, it is switched onto the automatic landing 
system. It is assumed that the aircraft is receiving the necessary position data from 
the ground-based airport landing system and receives control -surface positions from 
transducers. This is considered time zero with respect to the nominal path. 

(9) The control system nulls the state errors and brings the aircraft onto the nom- 
inal flight path before reaching the decision altitude. 

RESULTS AND DISCUSSION 
Flight Paths 

In this study, the control system was tested along seven different nominal flight 
paths which are described in table 2. The paths differed in initial altitude, initial heading 
angle, and in angle of descent. Flight paths 1, 2, and 3 are depicted in figure 3. 

Flight path 1 is the easiest path to fly because its initial altitude is the highest and 
there is more time to null out initial errors in the state variables. Conversely, flight 
path 3 is the most difficult since its initial altitude is 150 m below flight path 2 and 300 m 
below flight path 1. The aircraft makes a descending turn to the right of 360°, 270°, 
and 180° for flight paths 1, 2, and 3, respectively, at which point it is headed north toward 
the runway. At this position (roll out) on the path, the aircraft is at a 100 -m altitude and 
it rolls to a wing -level attitude. It then descends in a straight line toward the runway. 
The simulation terminated at initiation of flare where off-nominal errors were recorded. 

Flight paths 1A, 2 A, and 3 A had the same initial altitude and heading, respectively, 
as flight paths 1, 2, and 3. However, they had a descent angle of 6° rather than 3°. 

These flight paths require more control than the first three and less time is available to 
null off-nominal errors. The rate of descent and roll angle are about doubled and the 
turning radius and flight time are about halved. The rollout altitude was maintained at 
100 m. 

Flight path 3B is a variation of 3A. The initial altitude and rollout altitude were 
raised 100 m. This had the effect of moving the helix back from the runway and about 
doubling the time spent on the straight -in portion of the path. 
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Design Considerations 

The control task presented here differs from conventional landing -approach 
schemes in that errors are taken with respect to a moving point in state space rather 
than a fixed line. Thus, even if the aircraft was flying on the correct helical path with 
the correct attitude, the fact that it was late or early in time would indicate to the con- 
trol system that there were state errors; in this context, the proposed guidance scheme 
can be considered four-dimensional. At crowded airports, more precise control of time 
sequencing should prove advantageous. 

Because the feedback gain matrix K is computed automatically, the designer’s 
job is to choose the weighting matrices M, Q, and R of equation (4) and then to eval- 
uate the subsequent performance of the control system. In fact, this becomes an itera- 
tive procedure where the designer is trying to find an optimum set of weighting param- 
eters. This task is not trivial, but by using diagonal weighting matrices, a functional 
relationship between their elements and the subsequent performance can be established. 

The control system was required to null out initial errors in the states, and to 
offset the effects of steady crosswinds. It was also desired that it not be sensitive to 
normal changes in the atmosphere and uncertainties in the aircraft's aerodynamics. 

Wind gusts and sensor noise were not considered in this study. Large initial state 
errors tend to demand large rapid control motions. In order to avoid this, state errors 
are weighted less by reducing the elements of Q and control excursions are weighted 
more by increasing the elements of R.^ However, the performance index does not 
include control rates and the R matrix influences the time integral rather than the 
magnitude of control excursion squared. Therefore, it is not unusual to get large control 
excursions for relatively short time with consequent high control rates. If the Q and 
R matrices are well chosen, the control system will null out the initial errors slowly 
in a manner which will cause no passenger discomfort. As the aircraft approaches the 
touchdown point, the influence of the M matrix tightens up the control and acceptable 
values of final state errors can be achieved. 

This approach worked fine until the aircraft was subjected to a crosswind. It then 
became apparent that control in crosswinds was a much more stringent requirement and 
that a tighter control was necessary. To reconcile the two requirements (wide bandwidth 
control for crosswinds and narrow bandwidth control for large initial errors) two 
approaches were taken. First, high gains were used, but rate and position limits were 
placed on the controls. This led to undesirable oscillatory behavior, and in some cases, 
instability. The rate limit on angle of attack seemed to be particularly troublesome. 

The second approach was to make the elements of Q and R time varying. In this 

^The initial choice of the M, Q, and R matrices was made following the sugges- 
tion on page 149 of reference 3. 
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method, the feedback gains are low initially, while the state errors are large. While the 
large errors are being nulled, the feedback gains are gradually increased and the system 
is reasonably tight when it approaches touchdown. This second method works well except 
that roll rate became excessive when the aircraft rolled out of its turn. This maneuver 
occurs at a 100-m altitude when high feedback gains are desirable. This problem was 
circumvented by imposing a rate limit on roll. The values of Q, R, and M which 
were finally settled upon are: 
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where the value of the elements are given by 
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These weights were arrived at experimentally by using the iterative process described 
above. 
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There are several constraints imposed upon the design. To avoid passenger dis- 
comfort, angular rates should not exceed 10° per sec and by federal regulation, roll angle 
cannot exceed 30°. Engine thrust can be throttled between 12.45 and 124.5 kN (2800 
and 28 000 lb) which correspond to thrust-mass ratios of 0.3 and 3 m/sec^. Final errors 
in the states will be a function of how well the weighting matrices are chosen, the magni- 
tude of the initial state errors, wind velocity and direction, and of the particular flight path 
chosen. After the wings-level maneuver the aircraft is headed north at about 62 m/sec, 
and since the runways are normally much longer than necessary for this aircraft, an error 
in x of ±62 m (i.e., approximately 1 sec) is not unreasonable. However, errors in y 
of greater than ±10 m are unacceptable. The fact that the runway is normally only 50 m 
wide and the wing span is approximately 30 m makes an error of 10 m undesirable even 
though some correction will be made during flare. Errors in velocity V and flight -path 
angle y are best viewed in terms of sink rate z (z = V sin y). This error should be 
kept within ±1.0 m/sec while a negative error (slower descent) is preferred. Error in 
heading should be kept within 5°. An error of that size can be nulled by a decrab maneu- 
ver. The most troublesome errors in this study turned out to be in y during the 6° 
approaches. 


Representative Landing Approach 

Time histories from a representative flight are shown in figure 4. The nominal 
curves are shown as dotted lines and the simulated flight (path 2A) as a continuous curve. 
This was a 6° approach with nominal initial heading due east and initial altitude of 550 m. 
The six states and three controls are shown, along with an overview of the XY-plane. 

The aircraft started wings level with a heading of 80° (10® error). It was 100 m higher 
than nominal and had more than a 300-m error in both x and y. The velocity was low 
by about 20 percent. The aircraft was near stall, and consequently, the angle of attack 
and throttle setting were set high for trim. These errors which are extreme were cho- 
sen to demonstrate the control -system performance. The plane was subjected to a 
15 -knot constant east wind. Immediately, the aircraft rolled with some overshoot to the 
23° nominal bank angle. The angle of attack was decreased and the throttle was cut back. 
This made the aircraft simultaneously pitch over and increase speed while initiating its 
curved descending turn to the right. The rollout at a 100-m altitude is apparent in the 
curves of y, if/, and <p. The entire approach lasted about 82 sec and rollout occurred 
at around 70 sec. The simulation was terminated when the plane reached an altitude of 
around 15 m just above where the flare would be executed. Most of the motion is smooth 
but roll -angle rate definitely reaches its 10°/sec limit and the thrust-mass ratio stays 
on its lower limit of 0.3 (T = 12.5 kN (2800 lb)) for much of the run. The overshoot in 
roll angle (-10°) at rollout is undesirable, but is much less apparent in the 3° approaches. 
It could, most likely, be considerably reduced by making the nominal roll angle a contin- 


27 


uous function. However, this would complicate (perhaps unjustifiably) the entire control 
system. At decision height, the aircraft was within 1 sec of the nominal. Final errors 
for this flight were -71 m, 0.0 m, -0.1 m, 3.9 m/sec, -0.8°, 0.1°, and -0.6 m/sec for x, 
y, z, V, y, 'p, and z, respectively. 

Representative Feedback Gains 

The elements of the feedback gain matrix K for flight path 2A are shown in fig- 
ure 5. The magnitude of the gains, in general, increase with time due to the time varia- 
tion of the weighting matrices Q and R. Some of the gains have a pronounced increase 
in the last several seconds of flight due to the influence of the M matrix. (See 
figs. 5(c), (g), and (k).) During the straight-in portion of the approach, 10 of the 18 gains 
are zero or can be considered zero. Of course, the more of these elements that are 
zero or near zero, the easier it is to implement the control system. In order to deter- 
mine the influence of a particular gain, one assumes a large reasonable error in the 
state and then calculates the resulting change in the control. If the change is insignifi- 
cant, then the element can be considered zero. For example, an error of 3 m in z 
near the end of the flight would only change the thrust-mass ratio by 0.06. Whereas an 
error in x of 60 m would cause the ratio to change by 1.2. Since these are both large 
reasonable errors at the termination of flight, the element K 33 should be considered 
zero. Using this approach, it can be seen that after rollout, thrust -mass ratio is a func- 
tion only of x and V, roll angle is a function of only y and ip, and angle of attack is 
a function of x, V, z, and y. This partial decoupling is useful to the designer in the 
iterative process of choosing the weighting matrices. Applying the same reasoning, five 
elements of K can be considered zero throughout the entire flight. They are K 23 , 

K 24 , K 25 , K 35 , and K 30 . Each of the time histories of the remaining 13 gains would 

have to be stored in the on-board computer. Eleven points with linear interpolation 
could easily approximate any of these functions. Some of them could be fitted by low 
order polynomials in time which would require even less parameter storage. It appears 
that the entire gain matrix and nominal path could be stored in less than 200 words of 
storage. This number could probably be reduced but it is important to realize that any 
parameter change in the nominal-state trajectory will affect these gains. Each airport 
using this system would probably have several different nominal approaches and these 
would also probably differ among airports. To be practical then, the computer must 
have access to the correct set of parameters out of many possible alternative sets. Sev- 
eral possible ways of accomplishing this are 

(1) All necessary parameter sets are stored in an on-board mass storage device. 

(2) The particular set of parameters needed is calculated on the airborne computer 
prior to landing approach. 
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(3) The parameter set is transmitted from the airport to the on-board computer 
when the landing approach is assigned. 

(4) The feedback gains could possibly be expressed as analytic functions of time 
and the five primary parameters of the nominal flight path. 

The feedback gains of figure 5 were calculated using the weighting matrices M, 

Q, and R as given in equations (30) to (32). The gains were used in the landing 
approach shown in figure 4. This choice of weighting matrices and resulting gains 
appears adequate for this flight path and control task. 

System Performance 

The control system was tested on each of the seven landing approaches described 

in table 2. Results from the first six flight paths are given in tables 3 to 8 and from 

flight path 3B in table 9. Each of these seven tables represent 25 simulated landing 

approaches. On each approach, the aircraft was subjected to a constant wind of 15 knots 

magnitude from a direction if/ . The direction changed for each approach and it was 

w 

generated as a random variable with uniform distribution between -180° and +180°. 

Each approach was started with off-nominal errors in all state variables. These initial 
errors were generated as random variables. They were normally (Gaussian) distributed 
with zero mean. The standard deviations selected for the six state errors in x, y, z, 

V, y, and \p, respectively, were 100 m, 100 m, 30 m, 3 m/sec, 1°, and 3°. Each flight 
path had a distinct random sequence of sets of initial state errors. 

The reason for employing randomness in the control task was that the system 
showed that it could be tuned to a particular task. Originally, only 90° crosswinds were 
used because they were thought to be a worse case. Then it was discovered that winds 
from the aft quarter occasionally gave more trouble. Because this aircraft has a high 
lift-to-drag ratio, these winds were particularly difficult on the steep approaches. Also, 
certain combinations of initial errors were easier to correct than others. Therefore, a 
strong wind randomly directed was used and the aircraft was initialized with large ran- 
dom off-nominal errors. The controls were initially set such that the wings were level 
(0 = 0) and the aircraft was trimmed. 

Table 3 presents the data for the 25 landing approaches along flight path 1. The 
mean, standard deviation, and range of initial errors are given for each state. The mea- 
sured statistics of the initial state errors differ from their ideal values because they 
represent a subset of a very long pseudo random sequence which has the selected statis- 
tics. The initial off-nominal values of the controls are not really errors, but show a 
variation in initial trim values. Roll angle is a constant -12.2° since the approach always 
starts wings level rather than rolled over at 12.2°. The final off-nominal errors taken 
just above the flare altitude show how well the control system did its work. Statistics 
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are given on the six state variables and sink rate z. As can be seen, the control system 
worked very well on flight path 1. Poor performance usually makes itself apparent in 
final errors in y and z. For these 25 runs, the approaching aircraft was never more 
than 1 m from the centerline of the runway and the final error in altitude was never 
greater than 1.1 m. Since the nominal sink rate is 3.24 m/sec, the maximum error in 
final altitude corresponds to about one -third of a second in time. This same observation 
can be made with regard to x since nominal x is approximately 62 m/sec and the 
maximum final error in x for the 25 runs was 24.7 m. Table 3 also indicates how much 
control was needed to fly the aircraft. For each of the 25 approaches, the largest off- 
nominal control excursion is identified whether it be positive or negative. The initial 
transient of each run was not included since the initial off-nominal values of at least roll 
angle would always be extreme for the flight and it would mask the desired information. 

So, for example, the average extreme for angle of attack was 4.6° greater than the nomi- 
nal and the extreme varied over the 25 flights from 3.3° to 5.7° above nominal. The table 
also gives information on extreme control rates for each run. These are important 
mainly in how they affect passenger comfort. The statistics given are on the extremes 
of each flight including the initial transient. Finally, the table gives the number of flights 
which were unacceptable because one or more of the constraints, discussed in the section 
on design consideration, was violated. 

As can be seen from examining tables 3,4, and 5, the control system was quite 
effective on flight paths 1,2, and 3, respectively. The largest final errors for these 
75 approaches were 2.1 m in y and 1.2 m in z. There were no unacceptable flights; 
both control excursion and control rates were moderate. 

Tables 6 to 8 show the results for the 6° approaches. As can be seen, the control 
system did not do as well on these approaches. Final errors in y and 4 / were larger, 
and roll angle tended to have larger off-nominal excursions. 

On path 1A, the mean final error in y was a disturbing 9.3 m. Because the air- 
craft did not execute the wings -level maneuver fast enough, it tended to go off to the right 
(east), and then it would roll to the left in order to get back on the centerline. On these 
6° approaches, there were only about 13 sec from initiation of rollout to the conclusion of 
the simulation. Since the variation of the error is small (1.3), performance on this 
flight path could probably be improved by starting the rollout a few seconds early. This 
is, in fact, what pilots have been observed to do when flying this approach. On flight 
path 2A, performance is slightly worse and nine of these had final errors in y which 
were greater than 10 m and consequently were unacceptable. The variation in final 
y-error was again small, indicating that a slight adjustment in the flight path could possi- 
bly bring the errors within ±4.4 m. Finally, on flight path 3A, the control system proved 
to be completely inadequate with 18 of the 25 flights having final errors in y which were 
greater than 10 m. It is possible that another choice of weights, specifically 6yQ and 
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Syjyj, might improve performance on the 6° flight paths. However, the principal problem 
seems to be one of time. Flight time decreases monotonically from 212 sec on flight 
path 1 to around 60 sec on flight path 3A. It appears that the 6° approaches do not have 
enough time to null out initial errors and the adverse effect of crosswinds. In order to 
support this premise, flight path 3A was modified such that initial altitude was 500 m 
instead of 400 m, and rollout was accomplished at 200 m instead of 100 m. This new 
flight path is labeled 3B. Its parameters are given in table 2 and its performance is 
reported in table 9. In order to facilitate comparison of table 9 and table 8, the exact 
same sequences of initial errors and wind directions were used. The flight time for 3B 
was about 15 sec longer than 3A, and all of this was added to the straight-in portion of the 
flight. The marked improvement in performance of 3B over 3A is apparent when table 9 
is compared with table 8. On 3B, 8 approaches exceeded the allowable limit of a 10-m 
final error in y; whereas on 3A, 18 approaches exceeded it. Since this improvement 
seems to be caused by the increased flying time, it may follow that doubling the initial 
nominal altitude and rollout altitude on the 6° approaches would lead to performance com- 
parable to that experienced on the 3° approaches. 

Finally, the control system was tested for sensitivity to uncertainties in the aero- 
dynamics of the aircraft and for normal variation in the atmosphere. A particular landing 
approach was used. Initial errors and wind direction were not varied. Twenty -five sim- 
ulation runs were made, each with a different atmospheric density function and each with 
a different set of values for the aerodynamic parameters in table 1. The results of these 
flights are given in table 10. All six aircraft parameters were allowed to vary randomly 
from run to run. The values given in table 1 were used as the mean and 4 percent of 
those values were used as the standard deviation. The distribution was normal. This 
resulted in a random variation in these parameters of about ±10 percent. Air density, as 
described in equation (7), was also given random variation. The multiplying coefficient 
1.22 of that equation was allowed to vary normally with mean 1.22 and standard deviation 
of 0.0244. This resulted in a random variation of air density of about ±5 percent for any 
particular altitude. These variations seemed only to propagate into the final errors in 
x, z, and V and into the extreme control excursions in angle of attack. None of these 
final errors were large enough to cause concern. The control system appears to handle 
these normal variations and uncertainties quite well. 

CONCLUDING REMARKS 

Linear optimal regulator theory is applied to a nonlinear simulation of an aircraft 
performing a helical landing approach. The nonlinear equations of motion are developed 
and are linearized (time -varying coefficients) about a quasi -steady helical flight path. 

The nominal-state time histories are given as explicit functions of time and a numerical 
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method for determining the constant -control inputs is presented. Control of the system 
to the nominal -state trajectories is posed as a regulator problem with time -varying 
weighting matrices in the cost functional. A method of solving for the feedback gain 
matrix is reviewed. 

The theory was then implemented in a simulation of a small two-engined transport 
airplane, and system performance was measured for seven distinct approaches including 
flight path angles (descent angles) of 3° and 6°. On each approach, the aircraft was sub- 
jected to large errors in initial values of state variables and to strong steady crosswinds. 
The system was also tested for sensitivity to normal variations in atmosphere and to 
reasonable uncertainties in the parametric description of the aircraft. Statistical data on 
200 simulated landings are presented. 

The control system performed very well on all the 3° approaches and was reason- 
ably insensitive to changes in the atmosphere and to parametric changes in the aircraft 
model. Performance, which was measured in terms of terminal errors, violation of 
design constraints, and passenger comfort, was not nearly as good on the 6° approaches. 
These approaches are more difficult to execute since descent rate and roll angle are 
doubled, radius of the helix is halved, and flight time is halved. On the most difficult of 
these 6° approaches, the total flight time was about 60 sec, and performance was definitely 
unacceptable. Evidence is presented to support the contention that total flight time is the 
critical factor and that performance on the 6° paths could be greatly improved by starting 
the approach from a higher altitude which would result in a larger turning radius and 
longer flight times. This hypothesis needs to be tested further if 6° approaches are a 
requirement. 

The control -system design method used in this study is relatively straight forward 
and is easy to implement with the aid of a modern computer. The only difficulty is in 
choosing the weighting matrices for the cost functional. It was concluded that the weights 
should be time varying for the particular control task studied. This method should be 
extended to the more complex simulation, which includes both actuators and sensors and 
has six-degree-of -freedom dynamics, before being implemented on the aircraft. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
June 9, 1975 
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APPENDIX 

JACOBIAN MATRICES USED IN THE LINEARIZATION 
OF THE AIRCRAFT DYNAMICS 

The aircraft dynamics were linearized about a nominal-state trajectory and 
nominal -control time history. This was done in order to apply linear optimal regulator 
theory to the computation of the feedback gains. As is evident in equation (3) of the text, 
the Jacobian matrices 


A(t) = 91 
ax 


B(t) = |1 
v 7 au 


are needed. These are 
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-sin a cos 0 

mV 



mV 

V 

( — + T COS Q?) 
\da 1 

sin 0 

(^+ 

T sin a) cos $ 

sin a sin 4> 

mV cos y 



mV cos y 

V cos y 


(Al) 


where 


9L = qSCT - 5 — 

ao; 4 ^o/ iso 


IS - 2 « sc LHm) <“ - 
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APPENDIX 


and 


o o 
0 0 
0 0 

aF _ 0 0 

3X 

0 0 


0 

0 

0 

0 


cos y cos p 

cos y sin p 

sin y 

-2D 

mV 


0 (Ji^ sin a cos <t> - g cos 


y ~ m cos 


*r 


r-2 


! 

.0 0 0 


(L + T sin a ) sin p 
mV^ cos y 


-V sin y cos ip 
-V sin y sin \p 
V cos y 

g cos y 
£ 

~ sin y 

(L + T sin a ) sin <fr sin y 
mV cos ^ y 


-V cos y sin p 
V cos y cos p 
0 

0 

0 

0 


(A2) 


An assumption that atmospheric density is constant over the range of altitude considered 
is incorporated into these equations. This assumption is discussed in the text. It is not 
made in the nonlinear simulation which uses the feedback gains. 
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TABLE 1.- PARAMETERS OF THE AIRPLANE 
(GEAR DOWN, FLAPS AT 40°) 


Lift-curve slope Cl^, ra< ^”'*' 7.162 

Angle of attack for zero lift a 0 , deg -10.4 

Drag coefficient for zero lift Cp) j0 0.157 

Efficiency factor t] 0.0314 

Aeodynamic reference area S, 91,04 

Mass of vehicle m, Kg 40 823 



TABLE 2. - PARAMETRIC 


Flieht I^ial Rollout j Descent I Initial j 
JL . altitude altitude ; angle j heading i 
p n i ho, m -z SI , m J y, deg % deg ; 



CO 



Sink 

rate 

z, m/sec 

1 Radius 
: of helix, 
m 

\ Roll ! 
1 angle 
<f>, deg 

i 

! Turn 

1 . rate 
y>, deg/sec 

Rollout 
time 
tsi, sec 

Final 
time 
tf, sec 

3.24 

1822 

12.2 1 

1.95 

185.0 

215.5 

3.24 | 

1822 ! 

1 

12.2 ■ 

1.95 ; 

138.5 

| 169.5 

3.24 

j 1822 

12.2 | 

| 1.95 

92.5 

123.0 i 

6.48 

909 

23.3 ! 

3.89 

93.0 

! 108.0 

6.48 

i 

909 

23.3 

3.89 

69.5 

85.0 

6.48 

909 

23.3 

3.89 

46.5 

61.5 

6.48 

909 

23.3 

3.89 

46.5 

77.0 



TABLE 3.- STATISTICAL DATA FROM 25 LANDING APPROACHES 

ALONG FLIGHT PATH 1 


jy = 3°; z = 3.24 m/sec; Vw = 15 knots; p w = U(-180°, 180°); 
tf = 211.5 sec; hf = 13.7 m; Xf = -262 m| 


Variable 

Mean 

Standard deviation 

Range 




Initial off-nominal errors 



x, m 

-6.9 

64.7 

-146 

to 119 

y, m 

-13.0 

84.5 

-153 

to 193 

z, m 

-4.1 

32.1 

-85 

to 

76 

V, m/sec 

.6 

2.8 

-4.2 

to 

5.2 

y, deg 

.2 

.8 

-1.3 

to 

1.9 

P, deg 

.5 

3.1 

-6.1 

to 

8.6 

a, deg 

.1 

1.4 

-2.1 

to 

2.7 

<P, deg 

-12.2 

0 




T/m, m/sec^ 

-.05 

.14 

-.36 to 

.20 


Final off-nominal errors (at decision altitude) 


x, m 

-10.2 

9.5 

-24.7 

to 

-0.0 

y, rn 

.6 

.3 

.4 

to 

1.0 

z, m 

.0 

.6 

-.6 

to 

1.1 

V, m/sec 

.2 

.4 

-.3 

to 

.7 

7, deg 

.0 

.2 

-.3 

to 

.4 

P, deg 

-.2 

.1 

-.3 

to 

-.1 

z, m/sec 

.03 

.25 

-.28 to 

.40 


Extreme off-nominal excursions during each flight 

a, deg 

4.6 

0.5 

3.3 

to 

5.7 

P, deg 

-.4 

2.2 

-2.3 

to 

3.9 

T/m, m/sec^ 

.09 

.22 

-.40 to 

.40 


Extreme control rate during each flight 


a, deg/sec 

0.3 

1.5 

-2.5 

to 

3.3 

4>, deg/sec 

7.8 

3.1 

-6.7 

to 

9.4 

T/m, m/sec^ 

-.03 

.15 

-.23 to 

.25 


Number of unacceptable flights = 0 






TABLE 4.- STATISTICAL DATA FROM 25 LANDING APPROACHES 

ALONG FLIGHT PATH 2 


jy = 3°; z = 3.24 m/sec; V w = 15 knots; 0 W = U(-180°, 180°); 
tf = 165 sec; hf = 14.6 m; Xf = -279 mj 


Variable 

Mean 

Standard deviation 

Range 




Initial off-nominal errors 



x, m 

8.4 

118.1 

-292 

to 236 

y, rn 

10.5 

93.9 

-199 

to 179 

z, m 

2.3 

21.5 

-27 

to 

48 

V, m/sec 

.2 

2.7 

-4.8 

to 

6.9 

y, deg 

.2 

.9 

-1.4 

to 

2.9 

0, deg 

-.9 

2.2 

-5.7 

to 

3.0 

a , deg 

.1 

1.3 

-2.9 

to 

2.8 

0, deg 

-12.1 

0 




T/m, m/sec^ 

-.06 

.14 

-.46 to 

.19 


Final off-nominal errors (at decision altitude) 


x, m 

-9.4 

8.1 

-26.3 

to 

0.4 

y* m 

.8 

.4 

.3 

to 

1.3 

z, m 

-.1 

.6 

-.7 

to 

1.2 

V, m/sec 

.1 

.3 

-.3 

to 

.8 

y, deg 

-.0 

.2 

-.3 

to 

.4 

0, deg 

-.2 

.1 

-.3 

to 

-.1 

z, m/sec 

-.01 

.23 

-.33 to 

.41 


Extreme off-nominal excursions during each flight 

a, deg 

2.9 

3.2 

-3.7 

to 

6.5 

0, deg 

4.1 

.1 

4.0 

to 

4.3 

T/m, m/sec^ 

.08 

.19 

-.34 to 

.39 


E 

xtreme control rate during each flight 


a, deg/sec 

0.5 

1.9 

-3.0 

to 

4.8 

0, deg/sec 

8.8 

.7 

6.7 

to 

9.4 

T/m, m/sec^ 

-.0 

.15 

-.29 to 

.23 


Number of unacceptable flights = 0 
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TABLE 5. - STATISTICAL DATA FROM 25 LANDING APPROACHES 

ALONG FLIGHT PATH 3 

\y = 3°; z = 3.24 m/sec; V w = 15 knots; = U(-180°, 180°); 
tf = 118.5 sec; hf = 15.5 m; Xf = -296 m] 


Variable 

Mean 

Standard deviation 

Range 




Initial off-nominal errors 



x, m 

n.i 

95.2 

-176 

to 174 

y, m 

-12.8 

104.1 

-209 

to 217 

z, m 

-2.5 

30.8 

-48 

to 

73 

V, m/sec 

-.0 

2.1 

-3.7 

to 

5.2 

Y, deg 

-.1 

1.1 

-1.8 

to 

1.8 

P, deg 

.4 

2.5 

-4.1 

to 

6.2 

a, deg 

-.1 

1.0 

-2.4 

to 

1.8 

<P, deg 

-12.1 

0 




T/m, m/sec^ 

-.02 

..20 

-.35 to 

.29 


Final off-nominal errors (at decision altitude) 


x, m 

-9.4 

10.5 

-29.1 

to 

3.4 

y, m 

.8 

.6 

-.3 

to 

2.1 

z, m 

.1 

.5 

-.5 

to 

1.0 

V, m/sec 

.1 

.6 

-.7 

to 

1.2 

Y, deg 

.1 

.3 

-.3 

to 

.6 

ify, deg 

-.2 

.1 

-.5 

to 

-.1 

z, m/sec 

.09 

.21 

-.31 to 

.72 


Extreme off-nominal excursions during each flight 

a, deg 

2.8 

3.1 

-3.8 

to 

7.0 

<t> 9 deg 

7.1 

.2 

6.8 

to 

7.4 

T/m, m/sec^ 

.08 

.23 

-.18 to 

.49 


Extreme control rate during each flight 


a, deg/sec 

-0.2 

1.8 

-3.5 

to 

2.6 

0, deg/sec 

8.3 

.9 

6.1 

to 

9.4 

T/m, m/sec^ 

-.00 

.13 

-.24 to 

.19 


Number of unacceptable flights = 0 
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TABLE 6.- STATISTICAL DATA FROM 25 LANDING APPROACHES 

ALONG FLIGHT PATH 1A 

[y = 6°; z = 6.28 m/sec; V w = 15 knots; = U(-180°, 180°); 
tf = 105.5 sec; hf = 19.5 m; Xf = -186 m) 


Variable 

Mean 

Standard deviation 

Range 



Initial off-nominal errors 

x, m 

27.4 

112.7 

-219 to 276 

y> rn 

13.3 

106.0 

-224 to 237 

z, m 

-7.7 

26.7 

-60 to 40 

V, m/sec 

-.1 

3.6 

-7.2 to 6.4 

y, deg 

.1 

1.1 

-1.5 to 2.4 

ty, deg 

-.3 

3.3 

-6.8 to 6.3 

a, deg 

-.4 

1.9 

-3.4 to 3.7 

0, deg 

-23.2 

0 


T/m, m/sec^ 

-.07 

.11 

-.17 to .17 


Final off-nominal errors (at decision altitude) 

x, m 

-11.3 

11.3 

-29.3 to 12.1 

y> m 

9.3 

1.3 

6.6 to 11.6 

z, m 

-.5 

.6 

-1.1 to .8 

V, m/sec 

.1 

.6 

-1.4 to 1.1 

y, deg 

.1 

.2 

-.3 to .4 

V', deg 

-2.0 

.3 

-2.5 to -1.4 

z, m/sec 

.14 

.32 

-.47 to .62 


Extreme off-nominal excursions during each flight 

o', deg 

3.2 

3.8 

-3.7 to 8.4 

<P, deg 

19.1 

.2 

18.8 to 19.4 

T/m, m/sec^ 

.19 

.18 

-.17 to .47 


Extreme control rate during each flight 

a, deg/sec 

0.4 

2.1 

-3.3 to 4.1 

<p, deg/sec 

-10.0 

0 


T/m, m/sec 3 

.00 

.11 

-.21 to .22 


Number of unacceptable flights = 9 


41 



TABLE 7. - STATISTICAL DATA FROM 25 LANDING APPROACHES 

ALONG FLIGHT PATH 2A 

jy = 6°; z = 6.28 m/sec; V w = 15 knots; p w = U(-180°, 180°); 
tf = 82.5 sec; hf = 15.3 m; xf = -146 rnj 


Variable 

Mean 

Standard deviation 

Range 


Initial off-nominal errors 

x, m 

-9.7 

89.2 

-235 to 147 

y, rn 

20.8 

75.4 

-143 to 192 

z, m 

1.2 

26.2 

-41 to 56 

V, m/sec 

.3 

3.4 

-6.9 to 6.2 

Y, deg 

-.0 

.a 

' -2.1 to 1.4 

P, deg 

-.3 

2.5 

-3.5 to 6.9 

a, deg 

-.9 

1.7 

-3.6 to 3.2 

0, deg 

-23.2 

0 


T/m, m/sec^ 

-.07 

.11 

-.17 to .24 


Final off-nominal errors (at decision altitude) 

x, m 

-2.9 

10.3 

-29.4 to 15.8 

y> m 

9.4 

2.6 

5.3 to 14.0 

z, m 

-.4 

1.0 

-1.9 to 1.3 

V, m/sec 

-.5 

.6 

-1.5 to 1.1 

Y, deg 

.1 

.3 

-.4 to .4 

P, deg 

-1.8 

.5 

-2.8 to -1.0 

z, m/sec 

.02 

.24 

-.35 to .32 


Extreme off-nominal excursions during each flight 

a, deg 

3.3 

3.1 

-4.2 to 6.3 

<5 P , deg 

13.2 

.3 

12.8 to 14.0 

T/m, m/sec^ 

.08 

.16 

-.17 to .31 


Extreme control rate during each flight 

a, deg/sec 

-0.1 

2.0 

-3.4 to 3.8 

<p, deg/sec 

-10.0 

0 


T/m, m/sec^ 

-.03 

.11 

-.17 to .21 


Number of unacceptable flights = 9 
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TABLE 8.- STATISTICAL DATA FROM 25 LANDING APPROACHES 

ALONG FLIGHT PATH 3A 

jy = 6°; z = 6.28 m/sec; V w = 15 knots; = U(-180O, 180°); 
tf = 59.0 sec; hf = 15 m; xf = -140 m] 


Variable 

Mean 

Standard deviation 

Range 



Initial off-nominal errors 

x, m 

-35.6 

88.4 

-225 to 129 

y, m 

-22.2 

89.5 

-266 to 81 

z, m 

-5.2 

36.3 

-61 to 54 

V, m/sec 

.1 

2.5 

-6.0 to 4.2 

y, deg 

.1 

.9 

-1.3 to 2.3 

\p, deg 

-.3 

2.8 

-5.4 to 4.7 

a, deg 

-1.1 

1.3 

-3.0 to 2.4 

0, deg 

-23.2 

0 


T/m, m/sec^ 

-.07 

.09 

-.17 to .09 


Final off-nominal errors (at decision altitude) 

x, m 

-21.7 

17.9 

-57.4 to 11.4 

y, m 

13.2 

8.5 

-12.4 to 26.4 

z, m 

-.5 

.7 

-2.0 to 1.2 

V, m/sec 

1.1 

1.1 

-1.1 to 2.9 

r, deg 

.1 

.2 

-.4 to .3 

0, deg 

-2.5 

1.7 

-5.2 to 2.6 

z, m/sec 

.18 

.28 

-.34 to .63 


Extreme off-nominal excursions 

during each flight 

a, deg 

0.5 

4.2 

-4.2 to 7.2 

0, deg 

21.4 

.8 

19.6 to 22.8 

T/m, m/sec^ 

.23 

.20 

-.17 to .62 


Extreme control rate during each flight 

a, deg/sec 

i.i 

1.8 

-2.0 to 5.1 

0, deg/sec 

-10.0 

0 


T/m, m/sec^ 

.03 

.10 

-.15 to .15 


Number of unacceptable flights = 18 



TABLE 9.- STATISTICAL DATA FROM 25 LANDING APPROACHES 

ALONG FLIGHT PATH 3B 

\y = 60; z = 6.28 m/sec; V w = 15 knots; i// w = U(-180o, 1800); 
tf = 75.0 sec; hf = 14.7 m; Xf = -135 m] 


Variable 

Mean 

Standard deviation 

Range 



Initial off-nominal errors 

x ? m 

-35.6 

88.4 

-225 to 129 

y, m 

-22.2 

89.5 

-266 to 81 

z, m 

-5.2 

36.3 

-61 to 54 

V, m/sec 

.1 

2.5 

-6.0 to 4.2 

r, deg 

.1 

.9 

-1.3 to 2.3 

i//, deg 

-.3 

2.8 

-5.4 to 4.7 

, a, deg 

-.9 

1.3 

-2.9 to 2.6 

4>, deg 

-23.2 

0 


T/m, m/sec^ 

-.07 

.09 

-.17 to .09 


Final off-nominal errors (at decision altitude) 

x, m 

-17.1 

16.5 

-48.4 to 18.8 

y, in 

7.7 

4.0 

-3.1 to 14.2 

z, m 

-.6 

.6 

-1.5 to .7 

V, m/sec 

.8 

.8 

-.6 to 2.4 

y, deg 

.1 

.2 

-.7 to .4 

P, deg 

-1.2 

.6 

-2.2 to .5 

z, m/sec 

.14 

.3 

-.64 to .60 


Extreme off-nominal excursions during each flight 

a, deg 

0.9 

4.1 

-4.0 to 7.1 

<t>, deg 

21.3 

.7 

19.4 to 22.6 

T/m, m/sec^ 

.21 

.!9 

-.17 to .55 


Extreme control rate during each flight 

a, deg/sec 

1.0 

1.8 

-2.0 to 5.0 

4>, deg/sec 

-10.0 

0 


T/m, m/sec^ 

.0 

.1 

-.15 to .16 


Number of unacceptable flights = 8 



TABLE 10.- STATISTICAL DATA FROM 25 LANDING APPROACHES ALONG 
FLIGHT PATH 2 WITH VARIATIONS IN THE AERODYNAMIC 
PARAMETERS. RANDOM VARIATION OF 
APPROXIMATELY ±10 PERCENT IN 
AIRCRAFT PARAMETERS AND 
±5 PERCENT IN AIR DENSITY 


Jy = 3°; z = 3.24 m/sec; Vw = 15 knots; = 170°; 
tf = 165; hf = 14.6 m; xf = -279 m) 


Variable 

Mean J 

Standard deviation 

Range 



Initial off-nominal errors 

x, m 

186.0 



y, m 

157.9 



z, m 

-24.0 



V, m/sec 

1.4 



Y, deg 

-.2 



deg 

2.1 



a, deg 

-.6 



0, deg 

-12.1 



T/m, m/sec 2 

.0! 




Final off-nominal errors (at decision altitude) 

x, m 

1.7 

3.7 

-4.3 to 12.2 

y, m 

.4 

.1 

.3 to .5 

Z , 111 

1.0 

.4 

.0 to 1.6 

V, m/sec 

-.3 

.2 

-.8 to -.1 

r, deg 

.4 

.0 

.3 to .4 

deg 

-.1 

.0 

-.1 to -.1 

z, m/sec 

.37 

.05 

.27 to .46 


Extreme off-nominal excursions during each flight 

a, deg 

3.2 

2.6 

-3.8 to 5.1 

<f>, deg 

4.1 

.0 

4.1 to 4.2 

T/m, m/sec 2 

-.20 

.03 

-.34 to -.29 


Extreme control rate during each flight 

a, deg/sec 

-1.0 

0.1 

00 

o’ 

1 

0 

CM 

T— 1 

1 

cb, deg/sec 

8.9 

.0 

8.9 to 8.9 

T/m, m/sec 2 

-.18 

.00 

-.19 to -.18 


Number of unacceptable flights = 0 











Nominal values 
Simulated results 




Figure 4.- Concluded. 
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Figure 5. - Feedback gains for flight path 2A. y = 6 °; \pQ 


900; h() = 550 m . 
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Figure 5.- Continued. 
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Figure 5.- Continued. 
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Figure 5.- Continued. 
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Figure 5.- Concluded. 
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